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Abstract: 

It is known that orthogonal polynomials obey a 3 terms recursion relation, as well as 
a 2 x 2 differential system. Here, we give an explicit representation of the differential 
system in terms of the recursion coefficients. This result is a generalization of an 
expression of Fokas, Its and Kitaev for the symmetric case. 
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1 Introduction 



Orthogonal polynomials [jE^] continue to be extensively studied at the moment, partly 
because of their intrinsic interesting structures, and also because of their relationship 
with many physical and mathematical applications. For example they are the main 
tool for studying the large N asymptotics of the spectral statistics of random matrices 

0, d, i, il, in. 

Orthogonal polynomials are known to satisfy 3-term recurrence relations, as well 
as differential equations. These differential equations are very useful for the derivation 
of the large N asymptotics, and the associated Riemann-Hilbert problem 0. 

For any given weight of the form e~ y ( x \ where V(x) is a polynomial, finding the 
differential system is easy, but so far, no general expression was known (except when 
V(x) is an even polynomial |J). The calculations are straightforward but more and 
more involved as the degree of V increases. Here, we derive an explicit and concise 
expression, valid for any polynomial V(x). The expression is very similar to that of 

I- 

In addition, these orthogonal polynomials satisfy differential equations with respect 
to the coefficients of V(x), viewed as deformation parameters. Here, we derive an 
explicit expression for these deformation equations. 

Outline of the article: 

- In part 2, we introduce the orthogonal polynomials and set the notations. 

- In part 3, we give the expression of the differential system satisfied by the orthog- 
onal polynomials. 

- In part 4, we give the expression of the deformation equations satisfied by the 
orthogonal polynomials, and give an alternative proof of the results of part 2. 

2 Orthogonal polynomials 

Consider the family of monic polynomials p n (x) = x n + . . ., orthogonal with respect to 
the weight e~ v( - x ^: 

J^dxp n (x)p m (x) e~ v{x) = h n 5 nm , (1) 

where V(x) is referred to as the potential. The orthogonal polynomials exist under 
suitable asumptions on V, and are unique |J |12| . 

For simplicity, we will assume that V is a polynomial of even degree, with positive 
leading coefficient, and the the integration contour T is the real axisQ. 

2 Many of these asumptions can easily be lifted, see section 4. 
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Instead of the p n 's, it is more convenient to introduce the following quasi-polynomial 
functions: 



■Pn{x)e 



-hv(x) 



which are orthonormal: 



J dx 1p n (x)lp m (x) 



(2) 
(3) 



2.0.1 Recursion equation 



It is well known [Eg, |10[ that the orthogonal polynomials satisfy a three-term recursion 
relation, which reads: 



X1p n (x) = ^ Qn,m^m = 7n+l"0n+l + Mn + 7n0n-l , 



m=0 



where 



7n Qn,n- 



h n -i 



and f3 n — Q n ,r. 



(4) 



(5) 



The infinite matrix Q is symmetric {Q n ,m — Qm,n), and has only 3 bands {Qm,n = if 
\m — n| > 1): 

fPo 7i \ 

7i Pi 72 
72 Pi 73 

V- o ■-. ••. ••./ 



Q 



(6) 



2.0.2 The string equation 

The derivative of ip n can be expanded on the basis {ip m } m=0 ^ 

oo 
m=0 

Integrating Eq. (^) by parts implies that the infinite matrix P is antisymmetric: P n> . 
—P m ,ni and from p' n = np n _i + 0(x n ~ 2 ), one has: 



(7) 



Pn,m+2( V '(Q))n,m 
Pn,n-1 + \ (V'(Q)) n>n -i 



which implies 



P 



n,m 



-^{V'{Q)) n , 

+\ {V'{Q)) n ,, 



if m > n 

n 

In ' 



if m > n 
if m < n 



(9) 

(10) 
(11) 



This can be written: 

P = ~(V'(Q) + -V'(Q)-) , (12) 

where for any infinite matrix A, A + (resp. A_), denotes the upper (resp. lower) 
triangular part of A. 

In particular, Eq. (|lOl) implies the following equations, known as "string equations" : 
= V'(Q) n , n , — = V'(Q) n , n - 1 (13) 

2.1 Differential system 

The sum in Eq. (|7]) contains only a finite number of terms and can be written: 

1 degV" 

^'rSp) = E Vk,nipn+k - r) k , n -kipn-k where r) k , n = V'(Q) n , n+k . (14) 
1 k=i 

Note that 

degV 

V'(x)i) n (x)= ^2 Vk,n^n+k + Vk,n-k^n-k ■ (15) 
k=l 

Using Eq. (f|) recursively, one can write any with k > and any ipn-i-k with 
> as a linear combination of ip n and ^ n _i with coefficients polynomial in x, of 
degree k. Eq. ([14]) can thus be rewritten as a 2 x 2 differential system: 

dx v M%) J v M x ) J 

where T> n (x) is a 2 x 2 matrix with polynomial coefficients of degree at most d = deg V. 
It is known that it must have the form [H2~]: 



where m„(x) and v n (x) are polynomials of degrees degw n (x) < d — 2 and degt> n (x) < 
d — 1. What was not known so far, is how to express u n {x) and v n (x) in terms of the 
coefficients 7„ and (3 n (i.e. the coefficients of Q), apart form the case where V is an 
even polynomial ||. Here, we prove the following theorem: 

Theorem 2.1 The matrix T> n (x) is: 

I ( V'(Q)-V'(x) \ ( V'(Q)-V'(x) \ \ 

T, ( T \ -h/'( T \ f 1 °\+ ^ Q ~ x )n-l,n-l \ Q-* Jn-l,n\fO ~ln\ 

"n[x) - 2 v [x) y Q ( yi(Q)-yi{x) \ ( V(Q)-v'(x) \ J \j n J 

(18) 
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Proof: 

One can check from [d/dx, x] — 1 and from the compatibility with the shift (n — > n+1), 
that u n (x) and v n (x) must satisfy the following recurrences: 

= V'(x) + J n U n (x) + J n+ iU n+1 (x) + (/3 n - x)v n (x) (19) 

1 = {Pn - x)(j n+ iU n+ i(x) - J n U n (x)) + jl +1 V n+1 (x) - J^V n ^i(x) . (20) 

and these recurrences uniquely determine u n (x) and v n (x) for all n, provided that the 
initial terms are given. 

Now, consider the infinite matrix: 

Notice it is symmetric R n , m = R m ,n- 

Computing the diagonal (n,n) term of (Q — x)R(x) and using Eq. (|I~3"D gives: 



V'(Q)n, n - V\x) = -V\X) = ((3 n - x)R n>n {x) + 7n^n-l,n(x) + 7n+l#n+l,n (x) (22) 

and computing the (n, n ± 1) terms, using Eq. (O) gives: 



n + 1 

V r/ (Q)n,n+l = = (/?n - X)i2n n+ i(x) + 7„i? n _i in+ i(x) + 7 n+ i^ n+ i jn+ i(x) (23) 

7n+l 

Th 

V (Q) n ,n-1 = = (/?n " x)i? ni „_i (x) + 7 n i? n _i jn _i (x) + 7n+l-Rn+l,n-l (^) (24) 

In 

Using the symmetry, i? n _ l n+1 = i? n+l n _! can be eliminated: 

1 = (Pn ~ x)('J n+1 R nin+1 (x) - l n Rn,n-l{x)) + 'jl +1 R n+1 , n+ i(x) - 7^i? n _i in _i (x) (25) 

We have thus proven that u n (x) = R n ^ n -i(x) and v n (x) = R n ^ n (x) satisfy the recursion 
relations Eq. (|19|),Eq. (p0|). In order to complete the proof, one must verify that the 
expression is correct for n = 1. This is not difficult, and we don't write it here. 

In any case, we will give an alternative proof of this result in the next section, which 
does require the n = 1 case. 



3 Deformations with respect to the potential 

The deformation equations determine the variation of the orthogonal polynomials when 
the potential V(x) is varied. In this section, we don't assume anything about the 
potential (which need not be a polynomial), or about the integration contour (which 
can be any (non-connected) path in the complex plane; we don't have to integrate by 
parts). 
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For any integer k > 0, let the parameter Uk be such that: 



() V(x) = \x k . (26) 



diik k 
In particular if V(x) is a polynomial, we have: 

V {x) = £ -u k x k . (27) 
fc=i fc 

Expand the variation of ip n (x) with respect to Uk on the basis of {il> n }'- 

a oo 

^-M*) = E(^LiW ■ (28) 

a, "fc m=0 

The infinite matrix is antisymmetric (by differentiating Eq. (^)): 

{Uk)n,m = —{Uk)m,n ■ (29) 

It is easy to see that (take the derivative of Eq. (§) w.r.t u k , use dp n /duk = 0(x n_1 ) 
and use Eq. (^)): 

^ = -^((0*) + -(«*)_) • (3°) 
This implies that the infinite matrix Uk is finite band: (Uk) n ,m — if |n — m\ > k, and 
this means that Eq. fl2"E| ) is a finite sum. Notice that the diagonal part yields (from 
Eq. ©): 

Uk can be used to form a differential system (by expressing any ip m in terms of ip n and 
V> n _l with Eq. (§): 



where Uk{x) is a 2 x 2 matrix with polynomial coefficients of degree < k. For instance 
it is easy to compute IA\ [x) : 

U 1 (x) = --( (3n - 1 ~ X 2 \ ) . (33) 

Theorem 3.1 The expression for Uk{x) is: 

k k ( ( xh ~Q k \ ( xk -Q k \ 

7/ /_\ _ 1 ( x ~Qn-l,n-l A I | \ X ~Q )n,n-l \ X ~Q )n-l,n-. 

Uk[X) ~2k{ Q k njn -x k ) + k ln [ 

\ \ x Q / n, n V x Q / n,n~l 

(34) 

Notice that: 

d 

tiU k {x) = ——\n^ n (35) 
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Proof: 

Notation: for any infinite matrix A, let A\ denote the I th diagonal above (or below, 
if / < 0) the main diagonal. A + (resp. AJ) denotes the strictly upper (resp. strictly 
lower) triangular part of A. 

Since the matrix Q = Q~\ + Qo + Qi has only three non-vanishing diagonals, we 
have: 

(Q k+1 ) + = (QQ k ) + = (QiQ k ) + + (QoQ k ) + + (Q-iQ h ) + 

= Qi ((Q h ) + + {Q k ) ) + Q (Q k ) + + Q-i ((Q k ) + - (Q k ) 1 



= Q (Q k ) + + Qi (Q k ) - Q-i (Q k ) 1 • (36) 

This combined with a similar calculation for (Q k+1 ^j implies: 

gfc+l _ Q k+1 = Q{Qk+ _ Q k L) + Qi{Q k + Qkj _ Q _ i{Q k + ) (37) 



which, componentwise, reads (using Eq. (p0|)): 
d d 

- 2(k + 1)- 1pn(x) = ~2k- X1p n (x) + (Q k ) (7 n+ i^n+l ~ TnVW-l) 

OU k+1 3u k \ /n,n 



+ {ln(Q k ) -Tn+lfg") 

V \ / n,n— 1 \ / n,n+l/ 



Using Eq. 



9 9 
2(A;+1)- Vn(^) = -2k-—x^ n (x) 



du k+1 du k 

+ (Q k ) ((x-/? n )^-2 7n ^_i) 

\ / n,n 



.9 9 
-2(A; + 1)- Vn-i(z) = ~ 2k ^— 



+ (g fc ) (2 7n ^„ - (x - /?„_i)^„-i) 

V / n— l,n— 1 

+ f 7 n-l(g fc ) 1 -7n(Q fc ) (39) 

V \ ' n— l,n— 2 \ / n,n—lj 



which can be further simplified using 

+ lnQ n ,n-l ■ (40) 

We thus get a recursion formula for U k {x): 

—2{k + l)W fc+1 (x) + 2kxU k {x) = 
-(xQ k -Q k+1 ) \ f _ Q k Q k x 

V /n— l,n— 1 I n. / Vn-l,n Vn— l,n-l \ f/l-n 

o (*q* - Q k ^) ) +lln \ -Q k n, n Q^-i J ( } 

V / n,n ' 

from which the result follows easily. 
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3.1 Alternative proof of Theorem 

When V(x) is a polynomial, Eq. (O) implies that: 



and therefore 



deg V 

P = J2 Uk + ikU k (42) 

fe=0 



deg V 

T> n (x) = u k+1 kU k (x) . (43) 

k=l 



Using theorem 57TL as well as Eq. (|T3"D one immediately gets the proof of theorem [TI 



4 Conclusion and generalizations 

We have found a general expression of the linear differential system T> n (x), which is 
very convenient for explicit computations. 

The result was derived only for a polynomial V, with even degree and leading coef- 
ficient positive. It extends immediately to V of arbitrary degree and leading coefficient 
(but now the integration path is taken in the complex plane such that the integrals 
converge), provided that the orthogonal polynomials exist in that case. 

It could also be extended to a potential V whose derivative V is a rational func- 
tion. The only difference is that V'(Q) is no longer a finite band matrix, and computing 
((Q — o)~ 1 ) nn _ 1 or ((Q — a)~ l ) nn is not easy (other approaches like Geronimus trans- 
formations might be more efficient 0). 

Another generalization is when the integration path is not from — oo to oo, but has 
endpoints. Again, a differential system can be found, which has poles at the endpoints. 
The part without poles is still given by |2.1| , but the residues of the poles are not easily 
found with that method. Other approaches (such as the dressing method of Zakharov 
Shabat ]14|, [5]]) are more efficient. 



Another generalization is for biorthogonal polynomials. Theorems 2.1 and 3.1 



be generalized, and the expressions will be given in a forthcoming article fTJ . 
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